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Abstract 

We propose a nonlinear modification of the Schrodinger equation that possesses the main 
properties of this equation such as the Galilean invariance, the weak separability of composite 
systems, and the homogeneity in the wave function. The modification is derived from the rela- 
tivistic relation between the energy and momentum of free particle and, as such, it is the best 
relativistic extension of the Schrodinger equation that preserves the properties in question. The 
only change it effectively entails in the Schrodinger equation involves the conserved probability 
current. It is pointed out that it partially retains the linear superposition principle and that it 
can be used to model the process of decoherence. 
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Since its discovery in 1926, the Schrodinger equation, 



h 2 

has been subjected to a very thorough examination, both experimental and theoretical. In the very 
same year, Madelung showed that one can represent this equation as a system of two nonlinear 
equations, 



dt m 

^L A R --R(VS) 2 - 2hR^- - 2RV = 0, (3) 

m m v ' at 



for the amplitude R and the phase S of the wave function ^ = Re lS . They form the so-called 
hydrodynamic representation of the Schrodinger equation and can be obtained as the imaginary and 
real part of this equation, correspondingly. The first of these equations is the continuity equation with 
the probability current 

j Sch = -R 2 VS (4) 
m 

and the probability density p = R 2 . However, it is the other equation that is crucial for determining 
the stationary states of quantum-mechnical systems defined by the condition OR 2 /dt = 0. 

Let us recall that, heuristically, the free Schrodinger equation can be derived from the relativistic 
dispersion relation between the energy and momentum of a free particle 



E = c\jp 2 + (mc) 2 = mc 2 \j 1 + (-^-) ■ (5) 



mc/ 



One obtains it by discarding the rest energy term from the leading approximation to (5), 
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E = mc 2 + — , (6) 
2m' v ; 

and applying the first quantization procedure. To account for a wider range of physical situations, 
one can subsequently add the potential term to the RHS of (6). This leads directly to (1). 

It is the purpose of this paper to derive a relativistic extension of the Schrodinger equation which 
would stem from a closer approximation to (5) than the truncated first order approximation (6) and 
which, in fact, would be a nonlinear modification of the equation concerned.^ 

The Schrodinger equation is homogeneous of degree one in the wave function, however it is not 
necessarily clear if this property is underlied by some important physical principle or whether it 
entails any physically relevant consequences. It is thus conceivable that this property is an accidental 
feature of the Schrodinger equation and, thus, disposing of it in a nonlinear modification of this 
equation will most likely not cause any damage. Nevertheless, as pointed out in ||, nonhomogeneous 
modifications of the Schrodinger equation suffer from the ambiguity in the definition of the energy 
functional. Although this occurs also in homogeneous nonlinear variants of this equation, there exists 
a class of homogeneous modifications for which the quantum-mechanical energy functional defined as 

1 For a broader exposition of various aspects of nonlinear quantum mechanics see ||. More elaborate, though still 
far from comprehensive a list of approaches and motivations to modify this equation can be found in VM. 
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the expectation value of the corresponding Hamiltonian operator coincides with the energy defined 
as a conserved quantity within the Lagrangian framework ||. This indicates that the homogeneity 
is necessary but not sufficient to guarantee the uniqueness of the energy functional. Weinberg has 
elevated the discussed property to one of the fundamental assumptions of his original generalization 
of the linear framework of quantum mechanics [III fM. It was his hope, particularly expressed in 0, 
that the homogeneity of nonlinear variants of this equation might entail their separability^ Even if 
this is true in his modification and was shown to be valid in some others ||, [J, one can find examples 
of nonlinear variants of the Schrodinger equation that contradict this thesis ||. Yet, it should be 
noted || |7j that the homogeneity is essential for nonlinear modifications of this equation to possess a 
weakly separable multi-particle extension. 

The separability of noninteracting subsystems [§, [|] has much more profound and better understood 
physical implications. To cast more light on this issue, let us demonstrate the weak separability of the 
Schrodinger equation in the hydrodynamic formulation. We are considering a quantum system made 
up of two noninteracting subsystems in the sense that || 

V(x u x 2 ,t) = V 1 (x 1 ,t) + V 2 (x 2 ,t). (7) 

We will show that a solution of the Schrodinger equation for this system can be put in the form 
of the product of wave functions for individual subsystems for any t > 0, that is, ^(xi,x 2 ,t) = 
^i(xi, t)^/ 2 (x 2 , t) = Ri(xi, t)R 2 (x 2 , t)exp {i(Si(xi,t) + S 2 (x 2 , t))} and that this form entails the sepa- 
rability of the subsystems. The essential element here is that the subsystems are initially uncorrelated 
which is expressed by the fact that the total wave function is the product of t) and ^ 2 (x 2 , t) at 

t — 0. What we will show then is that the subsystems remain uncorrelated during the evolution and 
that, at the same time, they also remain separated. It is the additive form of the total potential that 
guarantees that no interaction between the subsystems occurs, ensuring that they remain uncorrelated 
during the evolution. However, such an interaction may, in principle, occur in nonlinear modifications 
of the Schrodinger equation even if the form of the potential itself does not imply that. This is due 
to a coupling that a nonlinear term usually causes between ^x(xx,t) and ^ 2 (x 2 ,t). As a result, even 
in the absence of forces the very existence of one of the particles affects the evolution of the other 
one, clearly violating causality. This particular feature of the Schrodinger equation distinguishes it 
from the Klein-Gordon equation. One arrives at the latter from (5) upon squaring it and employing 
the first quantization. The time derivative of this equation introduces a coupling between ^>i(xi,t) 
and ty 2 (x 2 ,t), rendering the equation inseparable in the sense discussed. Another physically signif- 
icant difference between these equations is in the probability density which can be negative for the 
Klein-Gordon equation. 

The discussed separability is called the weak separability since it assumes that the wave function 
of the total system is the product of the wave functions of its subsystems in contradistinction to the 
strong version of separability that does not employ this assumption. As shown by Liicke |J [TIJ , weakly 
separable modifications, such as the modification of Bialynicki-Birula || or the Doebner-Goldin mod- 
ification, QTT I can still violate separability when the compound wave function is not factorizable, and 



thus they are not strongly separable. An alternative effective approach to the strong separability has 
been proposed by Czachor HT2| . This approach treats the density matrix as the basic object subjected 



2 As stated in "The problem of dealing with separated systems has led other authors to limit possible non- 
linear terms in the Schrodinger equation to a logarithmic form;" and "The homogeneity assumption (2) makes this 
unnecessary." 
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to the quantum equations of motion which are modified^ compared to a nonlinear Schrodinger equa- 
tion for the pure state. It admits a large class of nonlinear modifications including those ruled out by 
the fundamentalist approach advocated by Liicke and even those that are not weakly separable as, 
for instance, the cubic nonlinear Schrodinger equation. 

The Schrodinger equation for the total system, assuming that the subsystems have the same mass 
m, reads now 



dR\R\ t h 
~~dt 



and 



+ 



711 



{(V, 



v 2 

?2 



R\R 2 2 (Vi£i + V 2 S 2 )]} 



dt 



dt 



URl^l + %Rl^ + ^Rlv, ■ (RI^SA + -R\V 2 ■ (R 2 2 V 2 S 2 

IJ- 111 V / 777 V 



h 



1 dR\ 
Rl dt 



m i<2 v 7 J J 



(81 



U - (A, + A 2 ) R X R 2 - 2h Rl R 2 d{Sl + S2) - * 
m ot m 

. „ h 2 



R X R 2 (Vi5i+V 2 5 ; 

7T7 V 



dS dS 

(V! + V 2 ) R X R 2 = —R^R, + —R X A 2 R 2 - 2hR 1 R 2 —± - 2kR 1 R 2 — 2 - 
m m at ot 

+—R X R 2 (VrSX + —R X R 2 (V 2 S 2 ) 2 - V 1 R 1 R 2 - V 2 R t R 2 = 



R\R 2 



- (V 2 S 2 ) 



ft 2 A1.R1 

2 

-Vo 



- 2h 



dSx | h 2 
dt m 



h 2 A 2 R 2 dS 2 t 

5 2h—+ 

m K 2 



dt 



0. 



(9) 



Implicit in the derivation of these equations is the fact that Vi/i • V 2 g 2 = 0, where f\ and g 2 are 
certain scalar functions defined on the configuration space of particle 1 and 2, correspondingly. What 
we have obtained is a system of two equations, each consisting of terms (in square brackets) that 
pertain to only one of the subsystems. By dividing the first equation by R 2 R 2 and the second one 
by R\R 2 , one completes the separation of the Schrodinger equation for the compound system into 
the equations for the subsystems. Moreover, we have also showed that indeed the product of wave 
functions of the subsystems evolves as the wave function of the total system. 

We would like the modified Schrodinger equation to possess the property of separability of compos- 
ite noninteracting systems in the sense described and to be Galilean invariant. It should also ensure a 
non-negative probability density emerging in the continuity equation. Moreover, we require that the 
probability current of this equation be conserved. These properties are certainly of physical relevance. 
As we will see, our approach will result in a homogeneous equation by way of construction, exactly as 
the Schrodinger equation is rendered homogenous when derived in the manner outlined above. 

3 What this means in practice is that the basic equation is the nonlinear von Neumann equation fl3|| for the density 
matrix instead of some nonlinear Schrodinger equation for the pure state. 
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To proceed, let us note that by expanding (5) in a series of p 2 and applying the first quantization 
we obtain terms \ 2n A n ^, where A c = — is the Compton wavelength. Employing the Schrodinger- 
Madelung representation, we observe that the terms concerned contain growing with n a number of 
terms involving derivatives of R and S. The essence of our approach is to extract from them the terms 
that would satisfy our criteria. It is easy to check that the following reccurence relation emerges for 
n > 1 



A n ^ 



(A n + iB n ) exp(iS) = { [(A - (VS) 2 )in-i - (AS + 2VS ■ V)5 n _i 
i [(A - (VS) 2 )B n ^ + (AS + 2VS • V)A,-i] } exp(iS-), 



(10) 



where A = R and B = 0. The first part of this equality is trivially valid for n = as well. 
It should be noted that by expanding (5) in a series of A n \l/ one obtains nonseparable terms, the 
only exception being the leading n — 1 term that gives rise to the Schrodinger equation. Moreover, 
such an expansion would not be Galilean invariant. To ensure that both of these properties are 
maintained in our modification, it is more convenient if not necessary to work in the Schrodinger- 
Madelung representation. To facilitate more general considerations, let us first examine in more detail 
the lowest order terms emerging from (10). These are 



A 1 = AR-(VS) 2 R = A-(VS) 



R, 



B l 



(AS + 2VS ■ V)R 



V • (R 2 VS 
R 



and 



A, 



A-(VS) 2 R- (AS + 2VS- V)- 



V- [R 2 VS 



A - (VSf 



v • (r 2 vs 

R 



+ 



R 



AS + 2VS ■ V 



A - (VSf 



R. 



(11) 
(12) 

(13) 
(14) 



Among them, only B\ and AB\ are separable, although not Galilean invariant. RB\ is proportional 
to the divergence of the probability current of the Schrodinger equation, but the continuity equation 
(2) is Galilean invariant only as a result of cancellation of noninvariant contributions coming from its 
both terms. To see this, let us recall that under the Galilean transformation, 



t — t', x = x' + vt, 

the amplitude of wave function does not change, the phase changes according to 

S(t, x) = S'(t f , x ) + mv ■ f + — v 2 t', 

2 

whereas the operators that affect these quantities in (2-3) transform as 

^ 9 d _ 



(15) 



(16) 



(17) 



As one can deduce from the general reccurence rule (10) and the last two equations, only the 
terms AA n and AB n could produce Galilean invariant and separable terms from the terms that 
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became separable in a preceding iteration, but were not necessarily Galilean invariant. The conditions 
of separability and Galilean invariance are stringent enough to select only one term that maintains 
separability in successive iterations of A and at the same time yields terms that are Galilean invariant. 
It is 

a = Bi = RAS + 2VS • VR. (18) 

Let us note that even if VS ■ Vi? is separable, it is no longer so when acted upon by RA as it gives 
rise to nonseparable terms 



RA(VS ■ VR) = 53 i?V • \diSV ■ (diR) + <9;i?V ■ faS) 



(19) 



The same applies to other terms of the form RA ^VA 2 ™^ • Vi?j (n > 1) that emerge in higher 
iterations. For this reason they will be discarded. Now, with the second term of (18) discarded, 
RAa = V • (ii^V 3 ^) + RARAS, but since RARAS is not separable, it should be excluded from our 
construction as should other terms of the form RARA n S (n > 1) appearing in subsequent iterations. 
Repeating this procedure for RA n a (n > 1) one observes that it effectively produces V • (R 2 VA n S). 

The multiplication of all terms of the expansion in question by R is required in order to arrive 
at the correct form of the continuity equation. For the term B\ this has already been pointed out. 
Namely, it is RB\ = Ra and not B\ itself that gives rise to the divergence of the probability current. 
In the same manner, the terms RA n a (n > 1) contribute to the divergence of a new probability 
current in the modified continuity equation. We see that the only change our construction of the 
relativistically extended Schrodinger equation brings about is in this current, or, what amounts to 
the same, in the continuity equation of the Schrodinger-Madelung system (2-3). One can write this 
equation as 

an oo 

" mc 2 A 2 5> t (A*A) n a, (20) 



fi- 



at 



n=0 



where c n are the coefficients of the expansion. Multiplying both sides of it by R one arrives at 



or 2 

dt 



2h 



m 



2k 



X 2 c n c n RA n a = — ]T A 2 "c n V • (R 2 VA n S) 



m 



n=0 



2h. 



m 



R V ^ \fc n A n S 

n=0 



(21) 



The last sum simplifies to G(\ 2 C A)S ) where 



G(x) 



vT 



x 2 - 1 _ _ sin 2 ^^) 



the second equality being a straightforward consequence of relations y/1 — x 2 
cosy — 1 = —2 sin 2 |. Eventually, the continuity equation reads 



OR 2 
dt 



+ V • j RM = 0, 



where 



2h 



Jrm = --i? 2 V (G(X 2 C A)S) 



(22) 

cos ^sin -1 x^j and 

(23) 
(24) 
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Alternatively, one can write the modified Schrodinger equation as 



3* h 2 
in—- = — 
at 2m 



-A--(V ■Jrm) + ~(v -jsch 
P v ' P v 



(25) 



which is equivalent to 
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2m 



A + -V 

P 



pV (2G(A*A) -1)5 



(26) 



In general, the only notable exception to this rule being the modification of Bialynicki-Birula and 
Mycielski, the nonlinear modifications of the Schrodinger equation do not have the classical limit in the 
sense of the Ehrenfest theorem. The discussed variant of the modification presented is not exempted 
from this either. The nonlinear terms it introduces lead to some corrections to the Ehrenfest relations. 
We will now work out these corrections. For a general observable A one finds that 



d_ 

dt 



(27) 



where the nonlinear contribution is due to H^l = Hr + iHi, Hr and Hi representing the real 
and imaginary part of Hnl, respectively. The brackets <> denote the mean value of the quantity 
embraced. Specifying A for the position and momentum operators one obtains the general form of 
the modified Ehrenfest relations 

d (r) = (p) + I u (28) 



m 



dt 



d 
dt 



(P) 



where 



h 



2m 



W 



dVrpHi 



J dVp (2H T VS - VH R 



(29) 



(30) 
(31) 



In the derivation of the last formula it was assumed that / dVV (pHj) = which for homogeneous 
modifications as the one in question seems to be well justified 0. For the extension concerned these 
integrals are found to be 



h 



h / dVrV ■ 



pV (2G(A^A) - l) S 



h 2 

-— I dVV ■ 

m 



pV {2G(\ 2 C &) - l) S] VS. 



(32) 



(33) 



The presented modification is the best relativistic approximation to the Schrodinger equation 
that shares with it the properties of separability for composed systems, the Galilean invariance, 
and the homogeneity in the wave function. Moreover, the probability density of equation (25) (or 
(26)) is a non-negative function unlike that of the Klein-Gordon equation that originates within the 
same framework. What distinguishes the modification discussed from other nonlinear variants of 
the Schrodinger equation |8|, [U], |] that encompass the same properties is that they make additional 
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assumptions or employ principles beyond those specified here. These, essentially, are necessary to 
keep the number of new terms at a reasonably manageable level. Such a problem does not arise 
here as, in principle, there is only one new term, which certainly adds to the attractiveness of this 
proposal. What is even more, our modification does not introduce any new constants. This makes 
it truly unique in this respect for all other nonlinear modifications of the Schrodinger equation do 
involve some new parameters. 

The modification in question supports stationary states of the linear Schrodinger equation for 
which S = —Et/h + const without affecting their energy E. Some other characteristic solutions to it 
include the plane wave, but, interestingly enough, also ordinary Gaussian wave packets and coherent 
states. Since both of the latter are the result of superposing more elementary wave functions, each of 
which is a solution to the nonlinear equation, this means that, similarly as the Schrodinger equation, 
our modification allows, although only partially, for the linear superposition principle. It is easy to 
understand how this comes about here. The third spatial derivative of the phase of a wave-packet, 
for simplicity in one dimension (% = 1), 

Tntx^ \ t 
S = — m arctan — , (34) 

vanishes as do its higher order derivatives. Since the current correction of (26) starts with V 3 S in the 
expansion, this means that as far as wave packets are concerned, the new term does not affect the 
Schrodinger equation at all. The same holds true for coherent states for which AS = 0. The nonlinear 
modifications endowed with the discussed property, that we chose to call the weak nonlinearity, are 
rare (see |5|, [T3J] for the only other known examples of this kind) . 



As noted, the modification proposed employes only the most general and universal assumptions. 
Because of that, it seems to be perfectly suited to describe phenomena of some universal nature. As 
pointed out in [15|], one of the most important phenomena in the quantum realm, the process of deco- 



herence responsible for rendering classical features of the ultimately quantum world, can be modeled 
by Hamiltonians involving imaginary terms. This is exactly the case of our modification. It introduces 
only such terms and it does so in a minimal model-independent manner: the imaginary terms that 
play the role of "optical" potentials emerge naturally, being an integral part of the modification, and 
do not have to be put by hand. An attractive consequence of this is that the process in question 
can be thought of as caused by relativistic nonlinear corrections. One can hope that the modification 
proposed will provide a good model to investigate this process in greater detail. 



Acknowledgments 

I would like to thank Professor Pawel O. Mazur for bringing to my attention the work of Professor 
A. Staruszkiewicz that started my interest in nonlinear modifications of the Schrodinger equation and 
Professor G. A. Goldin for his interest in this paper. A correspondence with Professor Wolfgang Liicke 
concerning the problem of separability in nonlinear quantum mechanics and a stimulating exchange 
of correspondence with Dr. Marek Czachor about many issues of nonlinear quantum mechanics are 
also gratefully acknowledged. This work was partially supported by the NSF grant No. 13020 F167 
and the ONR grant R&T No. 3124141. 



8 



References 

[1] S. Weinberg, Ann. Phys. (N. Y.) 194, 336 (1989). 
[2] S. Weinberg, Phys. Rev. Lett. 62, 485 (1989). 
[3] W. Puszkarz, Energy Ambiguity in Nonlinear Quantum Mechanics, preprint |quant-ph / 9802001 
[4] G. A. Goldin and G. Svetlichny, J. Math. Phys. 35, 3322, (1995). 



[5] W. Puszkarz, Higher Order Modification of the Schrddinger Equation, preprint |quant-ph/9710007 . 



[6] W. Puszkarz, Nonlinear Phase Modification of the Schrddinger Equation, preprint |quant- 



ph/9710010 



[7] W. Puszkarz, Extension of the Staruszkiewicz Modification of the Schrddinger Equation, preprint 
|quant-ph/9710009 . 



I. Bialynicki-Birula and J. Mycielski, Ann. Phys. (N. Y.) 100, 62 (1976). 
[9] W. Liicke, Gisin Nonlocality of the Doebner- Goldin 2-P article Equation, preprint |quant- 



| ph/9710"033| . 

[10] W. Liicke, Nonlocality in Nonlinear Quantum Mechanics, preprint |quant-ph/99040i6| . 

[11] H.-D. Doebner and G. A. Goldin, Phys. Lett. 162A, 397 (1992); J. Phys. A 27, 1771 (1994). 

[12] M. Czachor, Phys. Rev. A. 57, 4122 (1998). 

[13] M. Czachor, M. Kuna, S. B. Leble, and J. Naundts, Nonlinear von Neumann-Type Equations, 
preprint |quant-ph / 99041 10| . 

[14] A. Staruszkiewicz, Acta Phys. Pol. B14, 907 (1983). 

[15] M. B. Mensky, Phys. Lett. A 196, 159 (1995). 



9 



